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Abstract 

Due to the nonperturbative contribution to the hyperfine sphtting the mass 
of the n^Pi state is strongly correlated with the center of gravity Mcog{rfiPj) 
of the n^Pj multiplet: M(n^Pi) is less than Mco^{n^Pj) by about 40 MeV 
(20 MeV) for the IP (2P) state. For 6i(1235) the agreement with experiment 
is reached only if ao(980) belongs to the l^Pj multiplet. The predicted mass 
of 61(2^ Pi) is PS 1620 MeV. For the isoscalar meson a correlation between the 
mass of /ii(1170) (/ii(1380)) and Afcog(l^-Pj) composed from light (strange) 
quarks also takes place. 



I. INTRODUCTION 

Since the discovery of the he meson the hyperfine (HF) sphttings of the P-wave states 
in heavy quarkonia were investigated in many papers [0- ||^. In Refs. 0- [^] it was clarified 
why the HF shift of the he meson with respect to the center of gravity Mcogi^Pj) of the Xc 
mesons turns out to be small, Ahf(^c) = — 0.87±0.24 MeV |]^. It is due to a cancellation of 
the perturbative and nonperturbative contributions which are both small and have opposite 
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signs: Agp(cc) ^ -1.7 ± 0.3 MeV and Ag^(cc) ^ 1 MeV. Here the total HF shift Arf is 
defined in the following way 



For light mesons the HF splittings of the P-wave states are of special interest, since for 
them the perturbative spin-spin interaction is suppressed as for any L = 1 state, while the 
nonperturbative HF interaction is expected to become larger. In our study it will be shown 
that the nonperturbative contribution A^p, defined through the vacuum correlators, does 
dominate and Ahf(IP) is about 30 MeV. Although the magnitude of the splitting depends 
on such vacuum characteristics as the gluon condensate G2 and the gluonic correlation length 
Tg, the total Ahf(^-P) turns out to be positive in all cases considered. 

In our calculations of the HF splittings we shall follow the approach developed in Ref. § 
where the spin-dependent interaction is considered as a perturbation and averaging the spin 
factors in a meson Green's function is performed without the expansion in inverse powers 
of quark masses, used in the usual treatment Therefore the spin-spin potential from 
Ref. can be used for massless quarks and the HF splittings appear to be proportional to 
[/io(nL)]^^ where ^^{nL) is the effective dynamical mass of a light quark which is defined by 
the extremum of the Hamiltonian deduced from the QCD Lagrangian. It is essential that 
fio{nL) depends on the quantum numbers of the state considered and is not small; for the 
nP meson containing a light quark and antiquark, fiQ{lP) ^ 0.40 GeV and ixo{2P) k, 0.52 
GeV and ii^ilP) = 454 MeV, /io(2P) = 566 MeV for the nP ss states. 

For the isovector IP mesons (61 (1235) and the ground states of the aj mesons) the 
calculated Ahf(IP) is 39(19) MeV for two different vacuum gluonic correlation lengths: 
Tg = 0.3(0.2) fm, and with the use of the experimental mass of 61 (1235) we obtain that 



where the theoretical error comes from the uncertainty in the value of the gluonic length Tg. 
From this result an important consequence follows, namely, the number (|1.2| ) is compatible 
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Ahf = M,o^{n'Pj) - M{n' Pi). 



(1.1) 



Mcog(l^Pj, / = 1) = 1258 ± 10 MeV, 



(1.2) 



with the experimental masses of the aj mesons (n = 1) only if ao(980) (but not ao(1450)) 
belongs to the isovector l^Pj multiplet, i.e. ao(980) is a usual qq state. 

For the 61 (2P) meson the mass M(6i(2P)) ^ 1620 MeV is predicted. The situation with 
the isoscalar P-wave mesons [hi and fj) is also discussed and a correlation between the 
masses of /ii(1170) and M^^^{1^ Pj)=l2Ab MeV for /o(980), /i(1285), /2(1270), as well as 
between the mass of /ii(1380) and M^og{l^Pj) ^ 1420 MeV for /o(1370), /i(1420), /2(1430) 
(or Mcog = 1470 MeV if /2(1525) belongs to a multiplet composed of a strange quark and an- 
tiquark) can also be interpreted as a manifestation of a positive (~ 30 MeV) nonperturbative 
HF splitting. 



II. NONPERTURBATIVE HYPERFINE INTERACTION 



The HF splitting of the P-wave mesons originates both from perturbative and nonper- 
turbative interactions: 



AHF(nP) = A^^p (nP) + A^^inP) 



(2.1) 



where the perturbative term for L = 1 exists only in second order of as and will be discussed 
in Sec. 5. The quantity A^p is defined by the nonperturbative spin-spin potential which is 
usually presented in the form, 
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(2.2) 



As was shown in Ref. the spin-spin potential V^^{r) appears to be the same for heavy 
and light mesons (if the spin- dependent interaction is considered as a perturbation) and 
can be expressed through the vacuum correlators D{x) and -Di(x) which were introduced in 
Ref. ^ and calculated in lattice QCD Q- [|12| 



vrir) 



00 

2 J du 





3P'(r, u) + 3Di{r, u) + 2r 



(2.3) 



By definition, at the origin (x = 0) these correlators are related to the gluon condensate 
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D(0) + Di(0) = ^G2, (2.4) 

where the physical value of G2 = 0.04 ± 0.02 GeV^ is usually taken. 

In lattice calculations it was found that D{x) and Di{x) can be parametrized as expo- 



nentials at separations x > 0.2 fm: [|TT|- 113 



= (i exp ^-^j , = (ii exp l^-^^j , (x > 0.2 fm), (2.5) 

with the gluonic correlation lengths Tg and Tj^^ which turn out to be different in the quenched 
approximation and full QCD. In the general case the parameters d and d\, obtained in lattice 
measurements, differ from -D(O) and -Di(O). 

In full QCD with dynamical fermions (n/ = 4) the correlation length was found to be 
relatively large and the Di-correlator is small and can be neglected in some cases [p!2| : 

T, ^0.3fm, di^^d, (ny = 4). (2.6) 



It was shown in Ref. that in this case the correlator D{x) can be taken as an exponential 



over all distances, i.e. d = DiQ), 



X 



/^(x) = D(0) exp (^--J , (T.^O.Sfm) (2.7) 



and from Eq. (|2.4| ) in this case 



^(0) ~ T^^2 = 0.55G2. (2.8) 

lo 



(2.9) 



Then from Eq. ( |2.3| ) the potential V^^{r) is given by the expression: 

VfP(r) = 6rf Jexp(^-^^^^^^^rfz/ = 6rfriri(^^^, = D(0). 
The string tension a is defined in the general 

00 00 

a = 2 j dv j dAD(VA2 + z/2), (2.10) 


and for D{x) taken as an exponential at all distances it reduces to the relation 



a = .dT-, or rf=^, (2.11) 

If a is fixed and not large {a ^ 0.14 GeV^) then for the gluon condensate a reasonable value 
0.036 GeV^ (for Tg = 0.3 fm) follows. In this case the nonperturbative HF splitting is 



For light mesons the HF shift in the form of the relation (|2.12|) gives a dominant contribution 
also in cases when D{x) cannot be interpolated up to the origin, see below. The matrix 
elements in Eq. (|2.12|) will be calculated in our paper with the use of the solutions of the 
spinless Salpeter equation and the definition of the effective mass of a light quark will 
be discussed in the next section. 

Here we would like to notice that the potential V^^^(r) in Eq. ( |2.9| ), corresponding to the 



exponential correlator from Ref. [O, has an essential shortcoming. From our calculations 



it follows that this term gives a rather large nonperturbative shift in charmonium, 

A^^(1P, cc) > 5.0 MeV, {Tg = 0.3 fm), (2.13) 



so that the total splitting (|2.1j ) turns out to be positive for he in contradiction with the 
experimental negative number. Therefore, to explain the HF splitting of the IP state in 
charmonium one needs to know D{x) in detail at small distances, since the HF splitting in 
heavy quarkonia appears to be very sensitive to the behavior of the correlators D{x) and 
Di{x) at short distances (this problem will be considered in another paper). However, for 
the light P-wave mesons the behavior of the correlators D{x) and Di{x) at short distances 
was found to be inessential and for them the potential ^/'^(r) in the form of Eq. ( |2.9| ) can 
be used with 5-10% accuracy. 

Nevertheless, for completeness we give below expressions for the correlator D{x) and for 
Vl^^{r), modified such as to make clear that there exists the opportunity to combine a small, 
"physical" value of the gluonic condensate G2 and a small correlation length Tg. Otherwise 



the values fitted in lattice calculations (quenched approximation), Tg ^ 0.2 fm in Ref. |TT 



and Tg ^ 0.12 fm in Ref. |T3[, give rise to very large "unphysical" values of G2, ~ 0.14 GeV^ 
and 0.23 GeV^, respectively. 

To this end D{x) is supposed to be a constant at x < a;o which differs from from the 
coefficient d in Eq. (|2.5|) and can be taken as 



D{x) = constant = d exp 



Xo_ 



X < Xq, Xq ~ 0.2 fm, 



(2.14) 



while at X > D{x) is given by the exponential (|2.7|) as it was observed in lattice measure- 



ments. Then even for very small Tg = 0.6 GeV^ = 0.12 fm the small value G2 ~ 0.02 GeV 
can be obtained for the gluon condensate. For the modified correlator D{x), Eq. ( p. 141) , 
the modified nonperturbative spin-spin potential is 



vrir) 



6d 



e -'9\/Xi 



du exp 



T 



e{xo - r) 



+6drKi ( — j 9ir-xo). 



(2.15) 



For the P-wave light mesons the difference in the nonperturbative HF shift for the po- 



tential V^^{r) and V^^^(r) does not exceed 10% and therefore the simpler potential 



NP/ 



defined by Eq. (|2.9| ), can be used. Still for the he meson in charmonium such a modification 
of the spin-spin potential is important. 



III. SPECTRUM AND MATRIX ELEMENTS 

The fine structure and HF splittings in light mesons, with the exception of vr and K, 
are typically much smaller than their masses and therefore the spin-dependent interaction 
can be considered as a perturbation. Then the choice of an unperturbed Hamiltonian is of 
great importance and here the unperturbed approximation is formulated with the help of 
the spinless Salpeter equation, 

{2 v/p' + m2 + Voir)} ^„i(r) = E^l ^ndr), (3.1) 

where m is the current mass of a quark and Vo(r) is the static potential. We have chosen this 
equation since under some assumptions it can be deduced from the QCD Lagrangian. In 
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particular, if in the Feynmann-Schwinger representation |1T3|JI^ the backward trajectories 



are neglected, then for L = the QCD Hamiltonian for the spinless quark (antiquark) 
coincides with Eq. ( p.l| ) and for L = 1 the correction to the equation p.lj ) is not large [jl5[ . 
Therefore we can use the Salpeter equation for the P-wave states. 

For light mesons in Eq. ( |3.1| ) the current mass is taken to be zero and the static potential 
Vo(r) is taken in the form of the Cornell potential, 

Vo{r) = -^^ + ar + Co, (3.2) 
6 r 

where aes is an effective Coulomb constant. One can expect that for light mesons which 



have the rather large size -R > 1 fm, {R = y (r^)), the value of aes will probably be close to 
the so-called freezing value afj. = acs{r — > oo) which was found in Refs. [|l6l, [|l^, and has 
the value 

ttfr = 0.50 ± 0.05, (3.3) 

if the screening effects are neglected. However, even for such a large acfr, at long distances, 
r > 6 GeV~^, the Coulomb interaction is small compared to the linear confining potential 
and in most cases can be neglected. Therefore we consider here two variants: 

Qfeff = (caseA), aefr = 0.45 (caseB). (3.4) 



To fix the string tension a in the static potential (3^) one needs to take into account that 
although the Salpeter equation with a linear potential ar provides a linear Regge trajectory, 
however, as shown on Refs. [|I^], the slope of the Regge trajectory for the Salpeter equation 

a' = ^ (3.5) 
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differs from the slope a'^^ in the string picture where 

1 



(3.6) 



27r(Tst 

with the standard value of cxst ~ 0.182 GeV^. Therefore, to provide the experimentally 
observed slope the value of a in the Salpeter equation should be taken smaller than cXst: 



a = -(Tst = 0.143GeV^ (3.7) 

In most of our calculations just this number will be taken, but in some cases the value 
(T ~ (Jst ~ 0.18 GeV^ will be also used for comparison. Thus in case A the static interaction 
is characterized by the parameter a only, with its value given by the number (|3.7|) . With 
this smaller value of a the masses of the excited states in our calculations will be lower than 
in Ref. Ijl^ (where the same Salpeter equation was solved with Ust = 0.18 GeV^) and closer 
to the experimental meson masses for the excited states. 



IV. DYNAMICAL MASSES OF LIGHT QUARKS 

In Refs. IQ a relativistic Hamiltonian was derived from the meson Green's function in 
the Feynman-Schwinger representation with the use of the auxiliary field (einbein) approach. 
For L = and a spinless quark (antiquark) is given by the operator 

1 , ^ .1 



R 







where /^(r) and i/(/9) are the auxiliary operators and /x(r) is defined in the following way: 



.m4|. (4.2, 

In the definition ( [4.2|) r is the proper time and t is the actual time. With the use of the 
steepest descent method the extremal values /iex(''") = /^o and z/ex(/?) = i^o can be obtained 
with the following result: 



A^o = V + ^0 = (4-3) 

Then the relativistic Hamiltonian Hji in Eq. ( |4.1| ) reduces to the spinless Salpeter operator 

Hr = + /io + ar 2Jp + w? + ar (4.4) 

In what follows the extremal value /xo, which is an operator, will be replaced by the average 
of this operator which depends on the quantum numbers nL of the state considered, i.e. 
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fio{nL) = {^p"^ + m'^)nL form 7^0, 

fioinL) = {^)nL ifm = 0, (4.5) 

where m is the current mass of a quark (antiquark) and for hght quarks we take m = 0, 
while for the strange quark = 170 MeV will be used. 



The definition ([4.5|) of the effective mass of a light quark was already discussed in Ref. |T8 
where it was shown that the expectation value of Hf> in Eq. ( |4.4| ) coincides with that for 
the nonrelativistic Schrodinger Hamiltonian, if the effective mass is defined as in Eq. ([4.5|) . 

As seen from the definition ( |4.5| ) the dynamical mass of a light quark /io(nL) appears to 
coincide with half the average of the kinetic energy operator. 

/io(nL) = ^Ekin(nL) (4.6) 

In Table | the values of fio{nL) are given for different sets of the parameters of the static 
potential Vo(r). From Table | one can see that the influence of the Coulomb interaction is 
rather weak even for an a^s as large as a^s = 0.45, except for the IS" case, where it changes 
the dynamical mass by roughly 25%. This happens because the sizes of the light mesons are 
large, e.g. the root-mean-square radii R{nL) for the different states are as follows: 

i?(15) = 0.8^0.9fm; i?(25) = 1.3 ^ 1.4 fm; i?(35) = 1.6 ^ 1.8 fm; 
i?(45) = 1.9^2.1fm; i?(lP) = 1.0 ^ 1.2 fm; i?(2P) = 1.4 ^ 1.6 fm; 

i?(lD) = 1.3-^ 1.4fm; i?(2L)) = 1.6 ^ 1.8 fm (4.7) 

At such long distances the Coulomb interaction is small, only < 10% compared to the linear 
term ar. Moreover one can not exclude that at r > 1.2 fm screening of the Coulomb 
interaction may be important and therefore the Coulomb term in the static potential is even 
smaller and can be neglected, being important only for the 15* ground state. 

To illustrate our results, the spin-averaged masses of the low-lying mesons are presented 
in Table || and compared to the experimental values (isovector and isoscalar mesons) and 



TABLE I. The dynamical masses fiQ{nL) (in MeV) for different liglit mesons (the current mass 
m = 0). 





IS 


25 


35 


45 


IP 


2P 


ID 


Set A 

a=0.143 GeV^ 

Qeff = 


298 


445 


557 


650 


399 


516 


480 


Set B 

(7=0.143 GeV^ 

Oeff = 0.45 


375 


513 


616 


703 


436 


551 


508 


Set C 

(7=0.18 GeV2 

Qeflf = 


335 


500 


625 


729 


448 


579 


539 



also to the masses from the paper by Godfrey and Isgur [jT^, where the same Salpeter 
equation is solved for a different set of parameters: 



a = 0.18 GeV^ aGi(r) < a^r = 0.60, Co = -253, MeV m = 220 MeV. (4.8) 



As seen from Eq. ( |4.8|) in a rather large value was taken for the current mass m of a 



light quark, while in our calculations the best fit was obtained with Set A: 

a = 0.143GeV^ a^s = 0, m = 0, Cq = -357MeV. (4.9) 

The constant Co in Eq. (|]9D was chosen to fit Mcogi2^ Sj) = 1424 MeV. 

In Table |T| the experimental numbers refer to the isovector mesons which are not mixed 
with ss and are expected not to have a large hadronic shift. From this table one can see 
that 

(i) a better agreement with the experimental masses is obtained if ao(980) is a member 
of the l^Pj multiplet; 
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TABLE II. The spin-averaged masses Mcog(nL) (in MeV) of the low-lying light mesons. 



state 


2S 


35 


IP 


2P 


This paper 

a = 0.143 GeV^ 

"eff = 

Co = -357 MeV 


1424 
fit 


1870 


1241 


1707 


Ref. [17| 


1420 


1970 


1260 


1820 


experiment 
(/=1) 


1424 
±44 


> 1800 


1252") 
1306'') 


1632^=) 
1683^) 



") This value of Mcog(lP) is obtained if ao(980) belongs to the l^Pj multiplet. 

*) This value of Mcog(l-P) is obtained if ao(1450) belongs to the l^Pj multiplet. 

^'^ This value of Mcog(2^Pj) is obtained if 02(1660) belongs to the 2^Pj multiplet. 

'^^ This value of M(.og(2^P/) corresponds to the case when 02(1750) belongs to the 2^Pj 

multiplet. 
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TABLE III. The nonperturbative HF splittings Ag^(nL) (in MeV)for light mesons. 



state 


IS 


2S 


3S 


IP 


2P 


Tg = 0.3 fm 


125 


56 


30 


44 


27 


Tg = 0.2 fm 


96 


48 


25 


24 


20 



(ii) in our calculations the masses of the 35* and 2P states lie about 100 MeV lower than 
in |]T^ and are closer to the experimental numbers for Mcog(2aj) and 7r(1800). 

With the use of the dynamical masses fio{nL) = rriq, presented in Table |I|, the nonper- 
turbative HF splitting can be calculated, since from Eq. (|2.3| ) we obtain 



where we have taken into account the second correlator Di{x) in Eq. (4) to have the oppor- 
tunity to vary the values of the correlation length Tg. In particular for Tg = 0.2 fm the ratio 
di/d 1/3 was found in Ref . [ITT . 



In Eq. ( CT) 

J, = {rK,{r/Tg))^L, = {r K,{r /Tg)) - ^yK,{^)) (4.11) 

Here it is assumed that the gluonic correlation lengths Tg and Tj^^ in Eq. ( |2.3| ) are equal as 
it was observed in lattice measurements of D{x) and -Di(x) for = [|1^,|1^. We shall also 



fix the string tension a and from the definition Eq. ( p.lO| ) the parameter d is 

We estimate the accuracy of the calculated numbers to be about 10%. The nonperturbative 
HF splittings of the S-wave and P-wave light mesons are given in Table ^ for two values 
of the correlation length: Tg = 0.5 fm and Tg = 0.2 fm (in both cases a = 0.143 GeV^, 
aes = 0). 

As seen from Table |ITT| the nonperturbative HF shift is large, ~ 100 MeV, for the IS* 
ground state; for other states the numbers weakly depend on the value of Tg with the 
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exception of the IP state for which A^p is different for Tg ^ 0.3 fm and Tg ^ 0.2 fm, which 



are taken from the lattice measurements of the gluonic correlators |TT|- |T^. In most cases 
the magnitude of HF splitting is between 20^50 MeV. 

We consider also the P-wave mesons composed of a strange quark and antiquark taking 
for the current mass of a strange quark = 170 MeV. Then the dynamical mass of the s 
quark for different nL states turns out to be about 50 MeV higher than for a light quark (cf 
Table |) , in particular 

fio{2S,ss) = 505 MeV, /io(lP,ss) = 454 MeV, /io(2P,ss) = 566 MeV. (4.13) 

Correspondingly, the spin-averaged masses of the ss mesons appear to be about 170 MeV 



higher than those for light mesons, e.g. taking the set A of the parameters (|3.4| ) and the 
constant Cq = —250 MeV, defined from a fit to the spin-averaged mass of the 25* states 
(0(1680) and ?7(1440)), we have obtained that 

Mcog(lP, ss) = 1424 MeV, Mcog(2P, ss) = 1885 MeV. (4.14) 

At this point it is of interest to note that Mcog(lP, ss) coincides with the center of gravity of 
the multiplet: /o(1370), /i(1420), /2(1430) which are expected to have a large ss admixture, 
but it is 50 MeV smaller if /2(2P, ss) is identified with the /2(1525) meson. 

For the IP ss state the nonperturbative HF shift can be calculated from the expression 
i ^A% ) for = 0.3 fm and Eq. ( ^lOD for = 0.2 fm with the following result: 



AgP(lP,ss-) 



37 MeV, Tg = 0.3 fm. 



20 MeV, Tg = 0.2 fm. 



(4.15) 



V. PERTURBATIVE HYPERFINE SPLITTINGS 



From experiment it is known that the HF and fine structure splittings are practically 
small for all light mesons (with the exception of the vr and K mesons) small compared to 
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their masses and therefore the spin-dependent effects can be considered as a perturbation. 
Then, as was shown in Ref. the spin dependent potentials can be derived by averaging the 
spin factors which are present inside the meson Green's function defined in a gauge invariant 
way. In this approach the expansion in inverse quark masses is not used and in Ref. p[ it 
was deduced that to order all perturbative spin-dependent potentials Vi{r) {i = 1, 2, 3, 4) 
for light mesons coincide with those in heavy quarkonia with the only difference that the 
pole mass of a quark should be repaced by the dynamical mass jiQijiL) of a light quark 
(for a heavy quark fio{nL) coincides with the current mass to order as)- In particular, the 
perturbative spin-spin potential between a light quark and a light antiquark is defined as 

Then for the S'-wave mesons the perturbative HF splitting is given by the well-known ex- 
pression: 



where is the strong coupling in the MS renormalization scheme. In Ref. [T^ the spin- 

spin interaction was modified with a smearing function with a characteristic momentum 
scale of about 1.8 GeV. Consequently we can write in Eq. ( p.2|) for the 5- wave mesons 



~ as(l-8 MeV) = a,(M^) ^ 0.31 ^ 0.33. (5.3) 

Since the scale fi coincides with the mass M^- of the r-lepton we take here as{fi) = 0.31. 

The wave function at the origin entering Eq. (|5.2| ) cannot be precisely defined for the 
Salpeter equation, since the expansion of the wave function iprihif) (18) in a basis (which is 
used here for the numerical calculations as suggested in Ref. |T9|) is diverging at the point 
r = 0. Therefore we define Rno{0) = ip{nS, r = 0) as in the einbein approach p[ taking also 
into account the Coulomb interaction which gives a correction of about 10 -r- 20% and the 
largest one is for the ground state (^ 30%). Then Rno{0) can be presented in the form 



i?„o(0) = ^i2o{nS)a^{nS), (5.4) 
14 



where the coefficients C{nS) are the following: (ocfr = 0.39), ^(15*) = 1.31, ^(25') = 1.20, 
^(35) = 1.16, ^(45) = 1.14 and the values of the wave function at the origin are 

i?io(0) = 0.294 GeV^/^ i?2o(0) = 0.30 GeV^/^ 

i?3o(0) = 0.325 GeV^/^ i?4o(0) = 0.34 GeV^/^ (5.5) 

From these numbers one can see that the wave function at the origin is almost constant, but 
slowly growing because of the increase of the dynamical mass /io(nS') with n. 

The values of the perturbative splittings for the nS states are given in Table |I^ {c(jjs ~ 
as = 0.31). If one neglects the Coulomb correction in the wave function i?„o(0) then A^p 
will be about 30-50% smaller. To check our choice of -Rno(O) one can calculate the leptonic 
width of p(770): 

- _2«'|fi.«(0)IVi_« Y (5.6) 



Mj \ 37r 



which gives the following value for the leptonic width {o.^^ = 0.31; a = 1/137) 

re+e-(p(770)) = 7.36keV, (5.7) 

that turns out to be in good agreement with the experimental number re+e-(exp) = 6.77 ± 
0.32 keV [g (for = 0.33 the leptonic width is T^+e- = 6.8 keV). 

From the number (|5.5| ) for i?2o one can expect that re+e-(p(1450)) ~ 1.7 keV and the 
fraction T^+^~ /Ttotai for p(1450) is seven times smaller than for p(770). 

From the comparison of the nonperturbative and perturbative spin-spin splittings in 
Tables |III| and |V| one can see that for all nS'-states {n ^ 1) the perturbative splitting 
A^p(nS') turns out to be about two times larger than A^p while for the IS* state the 
nonperturbative contribution is larger. It is about 60% of A^p(lS'). 

Knowing the HF splittings we can calculate the masses of the isovector mesons (see 
Table [V]), neglecting the coupling to the other channels. 

We would like to notice here that all our calculations were done for a massless quark 
(antiquark) with only two parameters: the string tension a = 0.143 GeV^ (which defines the 
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TABLE IV. The hyperfine splittings of the S-wave Hght mesons (in MeV) with ajjg = 0.31 





15 


25 


35 


45 




194 


125 


94 


75(60) 


Ahf (total), Tg = 0.3 fm 


329 


185 


144 


96 


Ahf (total), Tg = 0.2 fm 


290 


173 


119 


95 


experiment 




165 ± 100 







TABLE V. The predicted masses of the 5-wave mesons in MeV [Tg = 0.2 fm). 





^(25) 


p(25) 


^(35) 


p(35) 


^(45) 


p{AS) 


theory 


1294 


1467 


1781 


1900 


2170 


2265 


Ref. 


1300 


1450 


1880 


2000 






experiment 


1300 ± 100 


1465±25 


1800±13 






2149±17 ''^ 



The mixing of 4^5*1 and 2^Di states is not taken into account. 
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dynamical mass of the quark (antiquark) /io(n5') and the spin- averaged spectrum) and the 
value (y.jjg ~ ajj^i^MT-) ^ 0.31 suggesting that the characteristic "smearing radius" is small 



as in Ref. |T^. Still, in such a simple picture the agreement with experiment is reasonably 
good and our masses for the 35* states are about 100 MeV lower than in Ref. [|r^ and close 
to the experimental mass of 7r(1800). 

To obtain the masses of the 45* states one needs to take into account the mixing of these 
states with the 2D states with Mcog(2D) = 1972 MeV (for the same set of parameters A). 
The mixing will be done elsewhere. 

VI. THE MASSES OF THE hi AND hi MESONS 

For the P-wave state the perturbative HF splitting is of order a1 and is expected to be 
small. To estimate the perturbative contribution one can use the expression ||20| 

„2 

3 ■K^ll{nP) 

This perturbative HF shift is negative and in ( |6.1| ) iriq is repaced by the dynamical mass of 
a light quark. This is allowed since the P-wave HF potential Vf{r) does neither depend on 
the renormalization scale nor on the mass of a quark (antiquark). This expression follows 
from the perturbative spin-spin potential for L 7^ 



9 7rm2 



9 a- — 

{r-')np ^ X ,"^„> {r-')np, (n, = 3). (6.1) 



y!ir)-laU\n,-'-W'^^^laU\-\n,] 1. (6.2) 



37r V3 ^ AJ r 37r V4 3 V r^' 

This short-range spin-spin potential has a characteristic size Phf which can be estimated 
from the value of the matrix element {r~^)np'- 

{r-''^)iP = 0.019 GeV^ {r-^)2P = 0.030 GeVl (6.3) 

If Phf(^^) = i{r-^)npy^^'^ then RuFilP) ~ 0.75 fm and Phf(2P) ^ 0.65 fm are rather 
large. From these estimates one can conclude that for the P-wave states Phf(^-P) ~ 0.65 fm 
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appears to be much larger than for the nS'-states where in the smearing function Ri{p{nS) = 
(1.8 GeV)^^ ~ 0.11 fm was taken from Ref. ||T^ At the distances i?HF ~ 0.65 fm the value of 



ajls needs to be taken at the smaller renormalization scale and is very close to the freezing 
value o.j^{q = 0) which is expected to be ajfsiq = 0) ~ 0.5. Therefore here we take 
ajjg{q = 0) ~ 0.45. The numbers obtained from Eq. ( |6.1| ) 

APp(lP) = -5.1 MeV, Ahf(2P) = -4.8 MeV, (6.4) 



are much smaller than the nonperturbative shift given in Table |ITl| and have opposite signs. 
Combining both contributions, one obtains the total HF splitting, 



Ahf(IP) 



39 MeV, ifT3 = 0.3fm 



19 MeV, ifTg = 0.2fm 



(6.5) 



or the average number Ahf = 29 ± 10 MeV. Knowing the mass of 6i(1235), 

M(6i(lP)) = 1229.5 ± 3.2 MeV, (6.6) 

the predicted mass for the center of gravity of the l^Pj multiplet {Tg = 0.3 fm) is 

Mcog(l^Pj) = 1258 ± 3.2 (exp) ± 10 (th) MeV. (6.7) 

The number obtained for Mcog(l^P7) is in surprisingly good agreement with the experimental 
mass Mcog(l^P7, exp) = 1252 MeV, if ao(980) belongs to the l^Pj multiplet, and does not 
agree with Mcog(l^Pj) = 1306 MeV obtained in the case that ao(1450) belongs to the l^Pj 
multiplet. Thus a strong correlation between the masses of Mcog(l'^Pj) and 6i(1235) follows 
from our analysis and to fit the experimental data one must assume that ao(980) belongs to 
the l^Pj multiplet and is a qq state. 

Then ao(1450) can be considered as a member of the 2^Pj multiplet with Mcog(2P) = 
1633 MeV from Table and therefore with the use of the total HF shift we predict for the 
mass of 6i(2P): 
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M{hi{2P)) = 1610 ^ 1618 MeV, (6.8) 

since the total HF shift from Table |Tl| and Eq. (|6^ ) is 

22MeVrg = 0.3 fm), 
Ahf(2P) = <i (6.9) 

ISMeVTg = 0.2 fm). 

In the approximation of closed channels used here the HF shift of /ii(1170) and 6i(1235) 
should be the same, see Eq. (|6.9|) . However, for /ii(1170) the experimental value of the HF 
shift is larger, 73 ± 19 MeV, and therefore one cannot exclude that /ii(1170) has a small 
hadronic shift, AMhad = 35 ±20 MeV (note that /ii(1170) has a much larger width, T{hi) ~ 
360 MeV, than 6i(1235)). There exists also the state /ii(1380) with M{^Pi) = 1386 ± 19 
MeV. It is assumed that /ii(1380) is mostly composed of a strange quark and antiquark. 
Then from the calculated Ahf (total) ^ 35 MeV {Tg = 0.3 fm and Agp = 4 MeV) one can 
obtain the center of gravity of the l^Pj multiplet of ss mesons: 

Mcog(l^Pj, ss) ^ M(l^Pi) + 35 MeV ^ 1425 ± 19 MeV. (6.10) 

This number can be compared with Mcog(l^Pj) obtained in the case if /o(1370), /i(1426), 
/2(1430) are members of the l^Pj multiplet and mostly ss states: 

M,(„^)(l^Pj) ^ 1422 MeV (6.11) 

and this experimental mass is in good agreement with the predicted mass ( |6.10|) . In the 
other case, when /2(1525) is a member of the l^Pj multiplet, the "experimental" value of 
the center of gravity, 

Mi^J^ifPj) ^ 1474 MeV (6.12) 

is not correlated with the mass of /ii(1380) and the shift of the mass of /ii(1380) appears to 
be larger (about 80 MeV) than in our calculations. 
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VII. CONCLUSIONS 



We investigated the nonperturbative spin-spin interaction in light mesons and estabhshed 
that 

1. For the IS* state the HF shift due to the nonperturbative effects is rather large, 
because the dynamical mass is relatively small, so that A^p ^ 0.4 Ahf(15', total), while for 
the excited nS states it is only about 15% of the total shift. 

2. Due to the positive sign of the nonperturbative HF splitting the mass of the n^Pi 
state is strongly correlated with Mcogin^Pj) being 30 ± 10 MeV smaller than Mcogin^Pj)- 
The value of this shift depends on the gluonic correlation length adopted. 

3. With the use of the mass of 6i(1235) our predicted mass of Mcog(l^-P7, / = 1) is 
1258 ± 10 MeV and this number is in agreement with the experimental masses of the aj{lP) 
mesons only if ao(980) belongs to the l^Pj multiplet. 

4. For 6i(2P) we predict the mass M(6i(2P)) ^ 1.62 GeV. 

5. Our analysis can be applied also to the isoscalar mesons where /ii(1170) and 
Mcog(l^Pj) = 1245 MeV lie rather close to each other if /o(980) is a member of the 1 Pj 
multiplet. 

6. In the approximation when /ii(1380), /o(1370), /i(1420), /2(1430) are considered to 
be composed mainly of a strange quark and antiquark, the difference A = M^ogi^^Pj, ss) 
— M(/ii(1380)) ~ 35 MeV is in full agreement with our estimate of the nonperturbative HF 
shift, A^p ^ 35 MeV for the correlation length Tg = 0.3 fm. 

7. The preferable value of the gluonic correlation length Tg = 0.3 fm was obtained from 
our analysis of the HF splittings of different mesons in accordance with the lattice data of 
Ref 

This paper was partly supported by the grant RFFI-00-02-17836. 
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